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qj ' Abstract 

. £h ! We employ the string/gauge theory correspondence to study properties of strongly coupled 

quark-gluon plasmas in thermal gauge theories with a large number of colors and flavors. In 
particular, we analyze non-critical string duals of conformal (S)QCD, as well as ten dimensional 
wrapped fivebrane duals of SQCD-like theories. We study general properties of the dual plas- 
mas, including the drag force exerted on a probe quark and the jet quenching parameter. We 
find that these plasma observables depend on the number of colors and flavors in the "QCD 
dual"; in particular, we find that the jet quenching parameter increases linearly with Nf/N c 
at leading order in the probe limit. In the ten dimensional case we find a non trivial drag 
coefficient but a vanishing jet quenching parameter. We comment on the relation of this result 
with total screening and argue that the same features are shared by all known plasmas dual 
to fivebranes in ten dimensions. We also construct new D5 black hole solutions with spherical 
horizon and show that they exhibit the same features. 
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1 Introduction 



The experimental data collected at RHIC [1], where heavy ions collide at very high 
energies, indicate that a new state of matter is formed during the collisions: the QCD 
quark gluon plasma (QGP). It partially behaves like a fluid with very low shear viscosity 
and, contrary to the common theoretical expectations, present evidence supports the 
picture of a strongly interacting medium [2]. Some of its observed features can be fitted 
with predictions from lattice gauge theory, but this cannot be used to describe its real 
time dynamical properties. 

A powerful tool to study strongly interacting gauge theories, both at finite and at zero 
temperature, is the string/gauge theory correspondence [3]. At present it has been used 
to study the properties of plasmas where gluons and adjoint matter fields (but no quarks) 
interact. The high temperature gauge theory at strong 't Hooft coupling g\ M N c 3> 1, 
and in the limit of large number of colors N c 3> 1, can be described by a dual classical 
supergravity model with background metric containing a black hole. The features of 
the metric depend on those of the dual field theory. The master example is provided 
by the thermal version of M = 4 Supersymmetric Yang-Mills (SYM) theory. The dual 
gravity background is AdS§ x S 5 with a large black hole in the center of the space. The 
above thermal field theory differs from thermal QCD in many aspects, one of the most 
evident being its lack of dynamical quarks. Nevertheless, the M = 4 SYM plasma is non 
supersymmetric and non confining as the QCD one, and it is tempting to compare the 
AdS/CFT predictions on J\f = 4 SYM and the extrapolated data from QCD in order to 
understand how macroscopic observables are sensitive to the microscopic details of these 
two theories. 

The string/gauge correspondence has been employed to evaluate the shear viscosity of 
the Af = 4 plasma [4] some years ago. This was a major success since the resulting value 
seems to obey a universal law and, furthermore, it is well compatible with the values for 
the QGP found at RHIC. This heavily supports the use of AdS/CFT to describe this 
physics. More recently, a lot of explicit calculations of real time dynamical properties, 
such as the drag force experienced by a quark moving through the plasma [5, 6, 7], the jet 
quenching parameter measuring the average energy loss by the parton per mean free path 
[8], the screening length of a moving color charge inside the plasma [9, 10], and the rate of 
photoproduction [11] were given. Generalizations of these results in different directions 
were pursued, respectively, in [12], [13], [14] and [15]. See, also, [16, 17, 18, 19, 20]. 

Despite some remarkable and unexpected numerical similarities between observables 
evaluated in the extrapolation to N c = 3 of large N c , M = 4 thermal SYM theory and 
the same observables deduced for the QCD quark-gluon plasma from the RHIC data, 
there are several reasons to be cautious. Of primary importance, it is mandatory to put 
the string/gauge theory correspondence at work for field theories closer and closer to 
real world QCD. In this sense, an unavoidable step to be able to properly talk about a 
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QGP consists in dealing with gauge theories coupled to Nf flavors of dynamical quarks. 
The first order effects in Nf/N c of the inclusion of flavors in the thermal string/gauge 
theory correspondence have been recently studied in [21, 22, 23, 24]. In this paper, we 
aim at considering the full, unquenched influence of dynamical quarks on the plasmas, 
by analyzing known supergravity solutions encoding these effects. 

In Section 2, we shall study plasmas having a non-critical string theory dual descrip- 
tion. The relevant gravity backgrounds we consider are AdS 5 and AdS 5 x S* 1 black hole 
solutions found in [25] . They are the non-extremal deformations of the zero temperature 
solutions found, respectively, in [26] and [27]. These are conjectured to be dual to confor- 
mal QCD and M = 1 SQCD in the Seiberg conformal window, in the so-called Veneziano 
limit [28], where Nf,N c — > oo but p = Nf/N c is kept finite. The gravity solutions are 
strongly curved and a' corrections are expected to be relevant. Nevertheless, it is worth 
exploring their properties, as they give a hint on the qualitative behavior of plasma ob- 
servables in relation with Nf. The holographic expressions for plasma properties like the 
screening length, the drag force coefficient, the jet quenching parameter and the diffusion 
coefficient for heavy quarks, can be immediately taken from the AdS critical cases, as 
they are not sensitive to the transverse directions. Interestingly enough, when expressed 
in terms of gauge theory quantities, these observables turn out to be independent on 
the number of colors and flavors in the "SQCD" model. This is not so for the "QCD" 
model, where, apart for the screening length which remains insensitive to the degrees 
of freedom of the plasma, there is an explicit dependence on N c , Nf. In particular, the 
drag force coefficient and the jet quenching parameter, which are the quantities on which 
we will mainly focus, are monotonic growing functions of p. This could be connected 
to a similar behavior of the entropy density. This qualitative behavior would reinstate 
the discussion about the use of the jet quenching parameter as a direct measure of the 
temperature of the plasma, since it would explicitly depend on the number of effective 
degrees of freedom. Finally, we make the (trivial) observation that the shear viscosity to 
entropy density ratio is still in these backgrounds, in the gravity approximation. 

In Section 3, we analyze a ten dimensional one-parameter family of black hole solutions 
found in [29] and conjectured to be dual to the high temperature regime of an M = 1 
SQCD-like theory with Nf = 2N C flavors and adjoint matter. The backgrounds describe 
the low energy dynamics of a stack of a large number of both wrapped "color" D5- 
branes and smeared "flavor" branes in the Veneziano limit. We study again the energy 
loss of a probe moving into the dual plasmas by evaluating the drag force and the jet 
quenching parameter. Moreover, we study the screening properties of the plasmas. We 
find that the drag force has analogous functional form as the one evaluated for the M = 4 
SYM plasma, that the jet quenching parameter vanishes and that probe quarks and anti- 
quarks are always screened. These results seem related to the critical thermodynamical 
properties of the model, in particular the fact that the black hole temperature, and so 
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the temperature of the gauge theory, coincides with the Hagedorn temperature of Little 
String Theory (LST). We argue that the same properties are common to all known 
models, with or without flavors (as the thermal Maldacena-Nunez solution [30, 31, 32]), 
having an UV LST completion. In particular, as shown explicitly in appendix B for the 
backgrounds found in [32], the LST plasmas are well expected to totally screen quark- 
antiquark interactions. The screening property of these plasmas should be responsible 
for the vanishing of the jet quenching parameter. 

We also present in appendix A a new family of black hole solutions which arise as 
compactifications on S* 3 of the Nf = 2N C configuration of [29]. Despite the inclusion of 
a new scale in the theory (the compactification radius), that naively might have been 
expected to modify the temperature away from the Hagedorn point, their holographic 
duals have the same thermodynamic problems and screening properties of the other LST 
plasmas. 

We conclude in Section 4 with a brief discussion. 

2 QGP and non-critical holography 

Non-critical string duals of 4d gauge theories with large N c , Nf both at zero and at 
high temperature have appeared in the literature in the last few years. In this sec- 
tion we investigate properties of the quark-gluon plasmas in these contexts by using 
the string/ gauge theory correspondence, though the gravity solutions at our disposal are 
generically strongly curved and a' corrections are not subleading. Our optimistic prej- 
udice, driven by the unexpected success of bottom-up approaches such as the so-called 
" AdS /QCD" duality [33, 34], is that the a' corrections may be not very large, and that 
the non-critical solutions should provide at least qualitative information on the dual field 
theories. 

The explicit backgrounds we will consider are an AdS§ black hole solution which was 
proposed as the non-critical string dual of the thermal version of conformal QCD and 
an AdS§ x S* 1 black hole solution conjectured to be dual to thermal M = 1 SQCD in 
the Seiberg conformal window [25]. The zero temperature models were constructed in 
[26] and [27] respectively. In both models the color degrees of freedom are introduced 
via N c D3-brane sources and the backreaction of Nf "flavor" branes on the background 
is taken into account. The flavor branes are spacetime filling brane-antibrane pairs, 1 so 
to reproduce the classical U(Nf) x U(Nf) flavor symmetry expected in the gauge duals 
with massless fundamental matter (see [36] for a very interesting analysis of this kind of 
setting). 

1 More precisely, in the six-dimensional model a brane-antibrane pair is just a single D5-brane doubly 
wrapped on the cigar [35] . 
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In order to study the dynamics of hard probes in the quark-gluon plasmas of these 
theories, we will make the assumption that the mass of the probes is related to the radial 
distance of a flavor brane from the center of the space, as it happens in the critical 
AdS$ x S 5 example, at least in some effective way. 2 In this case the general results on 
the drag force and jet quenching parameter obtained for AdSs black holes in the critical 
case, extend in a straightforward way to the non-critical setup. Indeed, the general form 
of the background string frame metric in the non-critical cases under study is 



ds 



2 _ fu\ 2 
RJ 



"// \ ,1,2 



U 4 



1 r dt + dxidxi 



dv 2 

+ (Ru) 2 - r ^ r + e^dQ k ; (2.1) 
■it — u H 



the dilaton = O is a constant. Let us now focus on the 5d model, postponing a brief 
analysis of the 6d one to the end of this section. 

2.1 QCD in the conformal window 

In the 5d model "dual to QCD in the conformal window", k = and the solution 
(in a' = 1 units) has the same radius, dilaton and five form field strength as the zero 
temperature one [26] 

R 2 = 200 

50 + 7p 2 - pv/200 + 49p 2 ' 



where 



= V200T4V-7P 

10Q C ' 1 ' 

F {5) = Q c Vo\(AdS) , (2.4) 

Q c and Qf are constants proportional to iV c and Nf respectively. 3 

Note that both the AdS radius (2.2) and the dilaton (2.3) - the latter being related 
to the gauge theory coupling by Qqcd ~ ~ depend on the flavor/color ratio p. Their 
dependence is the opposite: while R increases with p, the coupling decreases. The latter 
is consistent with the known fact that the zero temperature theory should be weakly 



2 Even in the model of [36], where also the branes supporting massive quarks are space-time filling, one 
expects a probe quark to be described by a string attached to the brane at finite distance, proportional 
to its mass, from the bottom of the space. 

3 If we use the same normalizations as in lOd, we have Qf = 2k§TiNf = iVy/(27r) 2 , and Q c = 
2k 5 T 3 N c — N c /(2tt). Even without the form of T 3 ,T4 but with only their ratio, one can get Qf/Q c = 
N f /{2nN c ). 
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coupled in the upper part of the conformal window, that is when p is the largest. Note 
that the behavior of the coupling is as expected in the Veneziano limit: it is given by 
jF(p)/iV c for some function T whose behavior is precisely J-{p) ~ 1/p for large p [28]. 
The black hole temperature and entropy density read 4 

uh_ A 3 ir 2 R 3 T 3 

where A 3 is the horizon area. The free energy can be obtained by suitably renormalizing 
the Euclidean action 5 



1 = w^Gi J d " x ^ e " 20 ( R + 4( ^ 0)2 + 5) " h F ^ ~ 2Qfe ~ 



(2.7) 



where 5 = 10 — d is the constant needed to cancel the central charge. Since the dilaton is 
constant, the DBI piece (the last one in (2.7)) plays the role of an additional cosmo logical 
constant term and the calculation is then practically the same as in ten dimensions [37]. 
Performing it, the free energy (density) reads 

The energy density can be derived directly as e = dpi ((3 = 1/T) or from the ADM 
definition of the energy 

E = J d 3 x^W\{K T - K ) , (2.9) 

where Kt, Kq are respectively the extrinsic curvatures of the 4d constant time black hole 
and zero temperature backgrounds. The result is the same: 

3n 2 R 3 T* , x 

e = ' (2 - 10) 

so that the thermodynamic relation 6 F = e — Ts is verified. The heat capacity (density) 
cy and the speed of sound v s are (see, for example, [24]) 

= Or, = , = - = 3 , (2.11) 

the latter result being consistent with the zero temperature theory being conformal. 



Here and in the following we use, for the d-dimensional Newton constant, the expression = 
(2^ 3 V) (d - 2)/2 .9 s 2 - 

5 We thank J.M. Pons and P. Talavcra for many discussions on this procedure and for common work 
on a related subject. 

6 And the relation s = —drF. 
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The holographic evaluation of the ratio between the shear viscosity r] and the entropy 
density in the plasma dual to this black hole gives the same universal result 

n - = f , (2.12) 
s An 

as for generic gauge theories equipped with gravity duals. The reason is that this black 
hole solution and the effective action from which it is deduced, have the general symme- 
tries advocated in [38] for this result to be valid. Formula (2.12) is somehow a trivial 
result in this setting, since we are limiting ourselves to the two derivative action, ignor- 
ing a 1 corrections that are most likely going to modify the ratio. If the latter were not 
the case, (2.12) would be a proof of the fact that the ubiquitous matter in the adjoint 
representation of critical string duals does not give any contribution to i]/s. Notice, in 
this respect, that the fact that matter in the fundamental representation does not change 
this result has already been argued in the ten dimensional model in [29] , that supposedly 
accounts for the full backreaction of the flavor branes in a SQCD-like theory. The result 
(2.12) is another example of this fact. From the expressions of the entropy density given 
above one can thus immediately deduce the shear viscosity 

ttR 3 T 3 , , 

All the thermodynamic quantities F,e,cv,r) have the same qualitative behavior as 
functions of p as the entropy density. The latter is a monotonically increasing function 
of Ql (for fixed p) and p. Its asymptotics in p are 7 

for p — > , 

(2.14) 

ffy/f (p 2 + CV)) forp^oc. 

The scalings with the number of colors Q c and the temperature T are the same as in 
M = 4, as expected. But, most importantly, this solution exhibits a nice dependence 
on the number of fundamental flavors Qf. As one increases Qf, s increases: the flavors 
contribute positively to the entropy density. At small p, (2.14), the correction to the 
pure glue scaling s ~ Q 2 C is the expected one, namely Q c Qf- Moreover, noting that the 
gauge theory 't Hooft coupling goes as A = e^°Q c ~ \[2 — 7p/10 + . . . for p — > 0, it is 
tempting to write the first p correction to the pure glue result as \Q c Qf. This is the 




7 Of course, these expressions must be taken with a grain of salt: since the background is expected 
to be corrected by order one terms, the numerical coefficients are not trustworthy. Moreover, since the 
dual field theory should be QCD in the conformal window for definite, finite values of p, the limit p — > 
(p — ► oo) is meaningful only as an indication of the behavior for small (large) p. The strict p = case 
is the finite temperature version of the unflavored Polyakov solution [39], that could be dual to a YM 
theory without flavors. 
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same enhancement effect of the flavor contribution found in the probe setting of M = 4 
SYM theory in [24]. On the other hand, the behavior for large p, (2.14), is found to be 
s ~ Q 2 and the first correction is already of order Q 2 . 

2.2 Drag force and jet quenching 

Now we turn to the discussion of the energy loss of partons in this plasma. As it is well 
known, the phenomenon of jet quenching observed at RHIC demands for a very efficient 
mechanism of energy loss. Recently, two main ways of accounting for this phenomenon 
have appeared in the literature. In [8], the main physical process for the energy loss 
relevant at RHIC, that is gluon radiation, is deduced in the perturbative QCD regime. 
All the effects of the strongly coupled plasma are encoded in just one quantity, the jet 
quenching parameter q, that accounts for the non-perturbative repeated soft interaction 
of the probe and the emitted gluons with the plasma. So, the problem is reduced to the 
calculation of q, which is postulated to be given by the coefficient of the L 2 term of a 
partially light-like Wilson loop of dimensions L~ ^> L. The calculation of this Wilson 
loop is then a straightforward exercise in string theory. This description should be the 
right one for partons with large transverse momentum. 

In a second approach initiated in [5, 6], the energy loss process is treated completely 
at strong coupling. The parton of mass m, moving with velocity v, is viewed as the 
extremum of a macroscopic string attached to a probe flavor brane. In order to keep 
the velocity fixed, one has to apply a drag force /, giving the parton some energy and 
momentum, that are ultimately lost in the plasma during the motion. This transfer of 
energy and momentum is seen as a flow of energy and momentum in the string, from 
the probe brane to the horizon. The configuration and the momentum flow of this string 
can be explicitly calculated, allowing to derive the drag coefficient p from the equation 
/ = pp, where p is the momentum. Then, of course, the coefficient enters the equation 
p = —pp, so its inverse is the relaxation time of the parton in the plasma. Being the 
entire process analyzed in string theory, and so at strong coupling, this description should 
be suitable for heavy quarks with relatively small transverse momentum. 

In this section, instead of calculating the two quantities p and q again (as we will 
do in the next section), we just note that, indeed, the calculations made for the critical 
AdS§ x S 5 black hole can be directly applied to the non-critical backgrounds. Thus, 
the general expression of the drag force (using the relation p/m = v/y/l — v 2 ) and jet 
quenching parameter are [5, 6, 8] 



dp 
~dt 




Vl^v 2 ' 



V 



(2.15) 




(2.16) 
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The novelty is that these quantities are now sensitive to the number of colors and flavors 
in our plasma. The asymptotics for the relaxation time are 



rn 




p + . . . for p — > , 



t ~ 



2nT 2 



*+°(*) 



(2.17) 



for p — > oo . 



Using as a representative temperature of the plasma T = 250 MeV, the relaxation time 
of a charm quark (taking m = 1.4 GeV) turns out to be 



t ~ 



0.7, 0.6, 0.5 fm 



for 



p = 0, 1, oo , 



(2.18) 



while t ~ 0.49, 0.40, 0.35 fm for T = 300 MeV and, respectively, the same values of p. 
Notice that the variation with p is very small and the values remain in the expected 
range of RHIC data (between approximately 0.3 and 8 fm [6]). 

The jet quenching parameter is a monotonically increasing function of p. Thus, it 
does depend on the number of degrees of freedom once the flavors are included. The 
asymptotics for the jet quenching parameter are inversely proportional to those of the 
relaxation time, 



If this qualitative behavior remains true in the full a' corrected background, this recasts 
the possibility of using q as a direct measure of the temperature of the quark-gluon 
plasma, as advocated in [8]. In fact, there is also a dependence on the "effective number 
of massless quarks" p; this is ultimately a function of the temperature too, albeit an 
unknown one. Putting in the numbers, 

q = 2.4, 2.9, 3.4 GeV 2 /fm for p = 0, 1, oo at T = 250 MeV . (2.20) 

In order to compare with [8], at T = 300 MeV one has 

q = 4.1, 5.0, 5.8 GeV 2 /fm for p = 0, 1, oo at T = 300 MeV , (2.21) 

so that it remains in the vicinity of the range of values compatible with RHIC data (5-15 
GeV 2 /fm [8]) for every value of p, and it is actually above 5.0 GeV 2 /fm for p > 1 and 
T = 300 MeV. Anyway, aside from the actual numbers, the interesting fact is that, as 
already noted, the variation of q is very small in the whole range of p. It signals the fact 
that the flavor contribution is not changing drastically the properties of the plasma. This 
is compatible with (and in a sense gives a reason for) the evidence that the actual values 




for p — > , 



for p — > oo . 



(2.19) 
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of plasma properties computed with gravity duals including only adjoint fields are very 
similar to the experimental ones. Even if at present the numerical values should not be 
taken very seriously as phenomenological predictions, their robustness suggests that the 
string/gauge correspondence provides an avenue that must be pursued. 

An interesting comment is in order at this point. It was shown in [20] that, for plasmas 
having lOd duals of the form AdS^-black hole times Sasaki-Einstein five-manifolds, ratios 
of jet quenching parameters are proportional to (the square root of) ratios of entropy 
densities. It was then conjectured that this result might also hold for approximately 
conformal systems such as the quark-gluon plasma of QCD. The authors of [20], however, 
ask themselves whether this relation is affected or not by the fact that some of the degrees 
of freedom in QCD are fundamentals. Let us assume that it is not affected. Then, it would 
be possible to extract the leading correction in Nj/N c to the jet quenching parameter first 
computed in [8] for M = 4 SYM theory, from the corresponding term in the expansion 
of the entropy density obtained for the D3/D7 system [23]. It gives 



a result that qualitatively agrees with the leading behavior (2.19) of our non-critical 
string/gauge dual prediction. Nevertheless, note that the conformal non-critical result, 
if not drastically modified by a' corrections, has a different behavior with respect to the 
proposal in [20] in that, say, the ratio of the jet quenching parameter to the one of M = 4 
SYM is not proportional to the square root of the ratio of the respective entropy densities 
at generic p. 

Let us point out that the holographic expression of other plasma properties like the 
screening length and the diffusion coefficient for heavy quarks, can be immediately taken 
from the critical AdS cases too, as any other quantity which is not affected by the 
transverse directions. Thus, in particular, from [9, 20], we conclude that the Debye 
screening length in our non critical setup will be only dependent on the temperature 
and the wind velocity, as in the M = 4 plasma case. Since at weak coupling the Debye 
screening length is known to depend both on the gauge coupling and on the number 
of colors and flavors, we are led to conclude, as in the M = 4 plasma case, that this 
could well be a genuine strong coupling effect. As for the above mentioned diffusion 
coefficient, it is related to the friction coefficient p by [5] D = T/pm; it is thus given by 
D = 2/(ttTR 2 ), and, then, it is a decreasing function of p in the present 5d setup. 

Let us conclude this part by warning the reader about the reliability of this non- 
critical string model. Aside from a' corrections, three other features have not been 
established so far [26]. The first one is a complete discussion of the stability: there 
are no explicitly tachyonic modes, apart from a possible closed string tachyon, whose 
behavior was not discussed. The second one is the exact spectrum of the D3-D4 system 
on a linear dilaton background in five dimensions (whose backreaction is supposed to give 




(2.22) 
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the AdS$ solution we have been discussing), that is not a priori guaranteed to contain 
(only) quark modes. Finally, if we take for granted the ten dimensional expressions for 
the mass of the open string tachyon (m 2 = —1/2) and the mass/conformal dimension 
relation (A = 2 ± y/A + m 2 R 2 ), the increasing of R 2 with p makes the dimension of 
the quark bilinear decrease in climbing the conformal window, contrary to the a priori 
expected behavior. 8 Of course, the assumption of the validity of the ten dimensional 
formulas is far from obvious. 

Despite this lack of control, this background seems to encode in a nice way many 
interesting features in a very simple setting, so we find it worthwhile to study them 
in a "phenomenologically" inspired way, guided by a similar spirit to that behind the 
" AdS /QCD" setting [33, 34]. In particular, this is the first case where the explicit (non- 
quenched) dependence on the number of flavors is exhibited for plasma observables in 
the holographic setting. 



2.3 SQCD in Seiberg's conformal window 

Let us briefly discuss the 6d model "dual to SQCD" (k — 1 in (2.1)), whose solution 
reads [27] 

R 2 = Q, ^° = A, e(t> ° = Jr> X = Qc6, (2.23) 

where x is the zero-form potential and 6 is the angle parameterizing Notice that the 
AdS radius R is now independent on both Q c , Qf, while the dilaton is not. As in the 5d 
model, the coupling has the correct behavior with p, but this time it is elementary. 
As in the 5d model, the temperature T and entropy density s can be computed, 

T = ^, s = 27n 2 Q 2 C pT 3 , (2.24) 

as well as the free energy F and energy densities e, 

F = J^Ql^\ e = 8i£!5kl!, (2 25) 

and the heat capacity cy, speed of sound v s and shear viscosity i], 

c v = 81n 2 Q 2 cP T\ vl = \, v= 27*QlpT* ^ 

In this case, the dependence of the thermodynamics quantities on Q 2 and p is just linear. 
Instead, both the drag force coefficient and the jet quenching parameters are independent 



8 We are assuming here that the open string tachyon is the gravity field dual to the quark bilinear, so 
that approaching the weakly coupled, higher edge of the conformal window its dimension should increase 
and tend to the free field value, namely three. 
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on Qc,Qf, since R 2 is a constant. The values are again slightly different from those in 
the M = 4 SYM theory. For a charm quark at T = 250 MeV, the relaxation time is 

0.5 fm, (2.27) 



irR 2 T 2 



that is smaller than the Af = 4 value but still in the allowed range that results from 
RHIC data. The jet quenching parameter at the same temperature is q = 3.6 GeV 2 /fm. 
In order to compare with the values in [8], at T = 300 MeV one finds q = 6.2 GeV 2 /fm, 
that is, higher than the M = 4 one and so more comfortably within the RHIC range. 

As in the 5d model, the screening length will be independent on the details of the 
plasma. The diffusion coefficient D will be independent on N c , Nf. 



3 QGP and wrapped fivebranes 

A one parameter family of black hole solutions in the background sourced by N c color 
wrapped D5-branes and Nf (smeared) flavor D5-branes was presented in [29]. This is 
conjectured to be the thermal deformation of the string/gauge dual to an M = 1 SQCD 
with quartic superpotential at the "conformal point" Nf = 2N C (f3 = 0, no breaking of 
U(1)r), coupled to Kaluza-Klein adjoint matter. The string frame metric and dilaton 
read (setting g s — 1) 



ds 2 = 



e^z 2 



Tdt 2 + dxl + N c a 



z 2 T 



dz 2 + - (dO 2 + sin 2 9 d^ 1 



1 ~ 1 



- — - (d9 2 + sin 2 9 dip 2 ) + J (dip + cos 9 dip + cos 9 dp) 2 



z 2 e^ 



T= 1 - ^ 



(3.28) 
(3.29) 



and £ is the parameter labelling the family. The horizon is at z ; for z = one recovers 
the zero temperature solution; t, Xi are dimensional, while the other coordinates are 
dimensionless. The temperature and entropy density of this BH are 



T 



1 



2ttVo 7 A 7 c 



AG 



(10) 



£(4-0 



(3.30) 



where A s is the area of the horizon. The fact that the temperature does not depend on 
Zq seems to be a common feature of black hole backgrounds of this kind, obtained from 
NS5 and D5-brane configurations. This also implies that T is independent on the energy 
density and, therefore, the free energy should be zero at this classical level. Then, the 
theory is in a Hagedorn phase and, indeed, the temperature above coincides with the 
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Hagedorn temperature of Little String Theory (LST). This suggests that there could be 
thermodynamical instabilities in this solution, in the very same way as it happens in the 
standard LST case [40, 41]. 

The instabilities are related to the UV LST completion of the model and thus cannot 
be cured by introducing IR cutoffs. In fact, we present in appendix A a new family of 
black hole solutions that should be string/gauge duals of the above plasma compactified 
on S* 3 , and argue that it should have analogous instabilities (as well as analogous screening 
properties) as its uncompactified counterpart. In view of these possible problems, it is 
not clear if these solutions can provide reliable descriptions of 4d finite temperature 
properties of string/gauge duals. Anyway, it is worth studying their features as long as 
they are the only known lOd black hole backgrounds where the backreaction of the flavor 
branes is taken into account. 

Concerning the hydrodynamical properties of the Nf = 2N C plasma, it was already 
pointed out in [29] that the ratio bewteen the shear viscosity and entropy density has the 
universal value 1/4-7T in the limit where the gravity approximation can be trusted. This 
result depends only on the symmetries [38] of the black hole solution and, thus, is not 
affected by the presence of the flavor DBI term in the gravity action (as it happens in 
[23]), nor by the troublesome thermodynamical behavior mentioned above. This seems 
to be related to the fact that the i]/s ratio depends on the properties of the horizon, 
while the LST behavior is an UV feature. 

3.1 The qq potential 

The goal of this section is to extract the potential between a probe quark and a probe 
antiquark in the Nf = 2N C QGP dual to the black hole backgrounds presented above. As 
in the zero temperature case studied in [29], we will find that the quark and antiquark 
are always screened. 

The Wilson loop for these two non-dynamical quarks separated by a distance L in 
the gauge theory coordinates should be evaluated [42, 43] by the action of a fundamental 
string that is parameterized by t = r, X\ — a, z = z(a). Using the string frame metric 
(3.28), we find that the length of this string is given by 

L(^)=4 V ^V> 7 f°° /V / ^~4 ^ = 2v /a^Ji-4k[4 

where is the minimum value of the radial coordinate reached by the string and Zq is 
the position of the horizon. By K[x] we denote the complete elliptic integral of the first 
kind. The length is a monotonic function that goes from a minimum value of 0, attained 
at z* = Zq, to the limiting value of it\/N~ c a', for — > oo. 



(3.31) 
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The energy of the string, upon subtraction of the infinite masses of the two quarks, 
reads: 



UK 
UK 



e 2 ^z 3 ^ 2 (z) 



dz — 



o *0 



z dz 



Z{) 



e*° zl 




1 - 



(3.32) 



It vanishes for z* — Zq «-> L — but is otherwise positive, which means that the quark 
and antiquark are always screened, as in the zero temperature case. As a matter of fact, 
we will argue in the following that total screening is a common feature of all known 
quark-gluon plasmas having LST as its UV completion. It is worth asking whether this 
affects or not other properties of the plasma, like the drag force on a probe quark and 
the jet quenching parameter. Let us address this question in the following subsections. 



3.2 Drag force 

Let us now evaluate the energy loss and the drag force of a heavy quark moving through 
the plasma, following the general procedure described in [5, 6, 16]. In the string/gauge 
picture, the quark is a macroscopic string with an extremum attached to a probe flavor 
D5-brane placed at distance z m from the black hole horizon. We consider a simple 
string configuration representing a test quark moving in a given spatial direction x: 
t = t, z = a, x = x(a, r). The static configuration x = constant corresponds to a quark 
at rest in the plasma. Its mass is given by the energy of the static string hanging down 
from z m to the horizon 

M rest = - — -/ y/-g U 9zz dz = — — -^(z^-zl). (3.33) 
2vra J Zo 2nVa' 

Thus, to get a very heavy quark we must require that the probe D5-brane is placed very 
far from the horizon. 

Let us now consider the string configuration corresponding to a quark which moves 
at constant velocity. This is achieved by setting x(a, t) = u(a) + vr. From the Nambu- 
Goto action on the background (3.28), it follows that the string Lagrangian density is 
proportional to 



z 1 



£ = e* x I AN C a' z 2 + z 4 T{z) u' 2 - AN C a' v 2 —-- , (3.34) 

J-\z) 
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from which it follows that the canonical momentum density conjugate to u(a) is 

». = t££m , A , (3.35, 

J4N C a! z 2 + ^ 4 T{z) u' 2 - 4iV c a' v 2 



Hz) 



and, finally, 



, V4N c a'A 2 \ z 2 (z4-z*-v*z*) 

For the above square-root to be real, at the point z c = zo/(l — v 2 ) 1 ^ 4 at which the 
numerator changes sign, also the denominator must change sign, this leading to 



e^ziv 2 



V 



2 



(3.37) 



The momentum loss of the steady moving quark is thus given by 

dp e* Zq v 

dt ~ 2na'^/T^v^ ' 

Assuming now that the following relation holds, at least for low velocity, 



(3.38) 



and that we can trust the relativistic relation 

M kin V 

P = R 1 , (3-40) 

where M kin is the kinetic mass of the quark, one obtains the coefficient 

»M kin = ^4 = ^e^z 2 N c T 2 , (3.41) 

or, equivalently, the relaxation time t = l//x. Let us notice that e* z 2 is the value of the 
dilaton at the horizon, e* • In the unflavored wrapped D5-brane setup of [30], the value 
of the dilaton at the minimal radial position e* multiplied by N c gives the parameter 
A = e* N c that measures the ratio of the string tension with the KK masses. This is the 
parameter that shall be fixed in the large N c limit, and must be much larger than one in 
order for the supergravity approximation to be reliable. Introducing also in our context 
9 the parameter A = e*° z 2 N c , we can write 

I^M km = 2n XT 2 . (3.42) 



3 Here A has not a clear interpretation as some cut-off coupling. 
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This is independent on the velocity and scales like T 2 , just as it happens in the M = 4 
SYM plasma [5, 6]. The only difference with the latter comes from the dependence on 
A, which is enhanced from A 1 / 2 to A (this is a usual phenomenon in passing from AdS to 
more general backgrounds with three-form RR fluxes). 

It thus seems that the energy loss of a probe quark moving in the QGP dual to (3.28) 
is not very much affected by the presence of dynamical massless quarks. 

Let us end this section by noticing that in order to have the friction coefficient /i 
disjointed from the kinematical mass Mfc in , we should analyze, as in [5], the small fluctua- 
tions around the trivial static configuration x = const., assuming then a time dependence 
of the type x ~ e~ M *. We do not pursue this analysis here. Instead we just notice that 
the limiting value of //, i.e., the one extrapolated in the extremal case where the probe 
brane touches the horizon z m — > zq, is expected to be as the one advocated in [16] for 
general asymptotically AdS backgrounds. This can be easily deduced from the general 
results in [19]. For a metric of the form 

ds 2 = -n(r) 2 A(r) 2 dt 2 + n(r) 2 [g^ dx^ dx v \ + fi(r) 2 [g„ dr 2 + gij dy i dy j ] , (3.43) 

where /i, v — 1, . . . , d; i, j — d + 2, . . . , D, the limiting value of \x is given by 



A** 



d r Vt 



2 \J cj TT Q — 



(3.44) 



r=r 



Applying this formula to the Nf = 2N C black hole, we find //* = ItxT as in the asymp- 
totically AdS cases. 10 



3.3 Jet quenching parameter 

Let us now evaluate the jet quenching parameter in the Nf = 2N C QGP under study, 
following the prescriptions given in [8]. What we need to compute is a rectangular light- 
like Wilson loop with one of the sides having a much smaller length, L, than that of 
the light-like side, L~. The prescription is that the Wilson loop in the fundamental 
representation (approximated by the renormalized Nambu-Goto action S^G,r) is related 
to the jet quenching parameter q, according to the relation 

Q = jpj2 SN G,r ■ (3-45) 

10 Incidentally, we notice that the same result is shared by all the other known (wrapped) fivebrane 
black hole solutions we will refer to in section 3.4. 
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Let us first rewrite the relevant part of the metric in light-cone coordinates x ± 
xJ/y/2: 



ds 2 = 



e*°z 2 



(F(z) + 1) dx+dx- + ^(1 - F(z)) {{dx + f + (dx') 2 ) 



+dxl + dx\ + N c a' \ 



z 2 T(z) 



dz 2 + 



(3.46) 



The string configuration for the loop is just x = r, x 2 = <J, z = z(cr), with the boundary- 
conditions z(±L/2) = oo. The Nambu-Goto action for a light-like length L~ reads, then, 



2L~e^z 2 f L/2 



S 

The equation of motion for x 2 then gives 

z' = 



Aa'N r z' 2 



1 + 



z 2 T(z) 



dz . 



(3.47) 



z 2 T{z) 



(3.48) 



where 7 is an integration constant, that gives 



4 v /]\r^ 7 



I = 



z dz 



zo 



(3.49) 



In order to evaluate the renormalized action, we have to subtract to this configuration the 
one corresponding to two separated strings extending from the boundary to the horizon: 



2L~ 
2n~a' 



POD 

J z 



9 zz = -^4r^V^-l ■ 



2^2n 



a' 



(3.50) 



We see that at small L, that is, large 7, we simply obtain S — S = 0, so that the jet 
quenching parameter q is zero 11 : 



C C C — 

*JNG,r — O — Dq 



2V2tc 



a: 



L I /2V4a' N c z' 2 I » r—— T 

-Ji+l ) -V^-i 



. (3.51) 



As we will point out in what follows, this is a general feature in gauge theory plasmas 
having Little String Theory as their UV completion. The q = result thus shares the 
same universality of the screening properties in LST plasmas. As such, these properties 
seem to be much more related to the 6d non local UV theories than to the IR local ones. 



11 This result was obtained in collaboration with C. Nunez. 
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3.4 On the general behavior of "LST plasmas" 

A possible puzzling result has emerged in studying the properties of the Nf = 2N C 
plasma obtained from a system of wrapped fivebranes: while the energy loss of a probe 
quark in the plasma, as described by the drag force calculation, has been found to be 
non zero, the jet quenching parameter turns out to be zero. It is known that the two 
quantities are related to the energy loss of a parton in two opposite regimes of transverse 
momentum (large momentum for the quenching parameter and small momentum for the 
drag coefficient), but obtaining such completely different results is an apparent puzzle. 
So, let us discuss this issue a bit. 

It is known that LST has a critical behavior at finite temperature: T = T H (Hage- 
dorn temperature) and negative specific heat characterize the BH solution related to 
unwrapped NS5-branes [40]. The same happens [41] to certain black hole backgrounds 
studied as finite temperature extensions of the wrapped D5-brane solution proposed in 
[30] as a dual to M = 1 SYM theory (the MN solution). It is not difficult to realize 
that both in the flat LST case and in the wrapped setups, the drag force is non zero, 
while q = 0. This also happens in the thermalized version [44] of the 3d setup studied in 
[45, 46]. 

Moreover there are more general BH versions of the unflavored MN solution, which 
were studied in [32]. There, a family of finite energy BH was found, whose temperature 
is above the Hagedorn temperature of LST and that basically have the same problems. 
Again, zero jet quenching and non zero drag force are expected to characterize these 
solutions. Checking the occurrence of total screening is not an immediate task, mainly 
due to the lack of analytic control on the solutions in [32]. As we show in appendix B, 
which can be taken as a reference for analogous calculations in general LST cases, these 
solutions manifest total screening too. As it transpires from the calculations in appendix 
B, it is expected that the same should happen for possible Nf ^ 2N C plasmas in the 
framework of [29], even though the correponding gravity backgrounds are not currently 
known. A possible way to evade the above critical behavior could have been to put the 
black hole on S* 3 . In this manner, another scale is introduced by hand into the system, 
and it might be expected that this allows to tune the temperature of the system away 
from the Hagedorn point. As we show in appendix A, this is actually not the case and 
these solutions behave in the same way as the other "LST backgrounds". 

These are indications that our "jet quenching puzzle" is not related to the presence 
of dynamical quarks in the lOd backgrounds. 12 Instead, let us observe that the jet 

12 There are phenomenological AdS/QCD models like [18] where dynamical quarks are argued to be 
present and responsible for the breaking of the flux tube connecting two probe quarks beyond some 
critical distance. The drag force and jet quenching parameter for those models were found to be both 
different from zero [17] with the jet quenching parameter going like T 4 . Together with the results of 
Section 2, this indicates that dynamical quarks seem to create no specific problem to the evaluation of 
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quenching parameter is definitely dependent on the UV behavior of the dual backgrounds, 
and so on the 6d LST asymptotics, while the same dependence for the drag force is not 
stringent (indeed we could consider a probe flavor brane "not so far" from the horizon, 
so to cut off the UV zone). 

This fact suggests that the result q = has to be associated more to non-local LST 
modes than to their claimed local counterparts. Analogous conclusions can be reached 
for the screening properties of the wrapped fivebrane plasmas, which again appear to be 
universal. In fact, the total screening result is also sensitive to the UV physics, and so 
to the LST regime. In conclusion, it seems reasonable to argue that quantities crucially 
sensitive to the UV behavior of the metric are determined by LST modes. These give total 
screening and, as a consequence, zero jet quenching parameter. Thus, these backgrounds 
do not seem to be useful in order to study UV properties of realistic plasmas. 

On the other hand, there are properties that seem not to be affected by almost any 
feature of the plasmas. For example, the r]/s ratio has the usual form because it is 
determined by the universal properties of the horizon. Analogously, the result of the 
drag force calculations indicates that the form of pM^ in is exactly the same in all the 
LST plasmas and practically the same as in M = 4. It is insensitive to the IR (and 
UV) details of the dual theory 13 and it is not substantially altered by the presence of 
flavors. In this respect, it seems that the indications it gives for the energy loss process 
in the dual field theory are very general and universal. The friction that a heavy parton 
experiences in a strongly coupled plasma is practically always the same, irrespectively of 
the features of the dual field theory. 

4 Concluding comments 

In this paper we have analyzed some properties of the known backgrounds dual to finite 
temperature 4d gauge theories with a large number of flavors Nf and colors N c , in the 
limit where their ratio p = Nf/N c is kept constant and different from zero. Our aim was 
to obtain information on the dependence of plasma features on the number of flavors, in 
order to better understand the connection of the string theory calculations with RHIC 
physics. 

In particular, we have studied plasma features of non-critical gravity backgrounds 
conjectured to be dual to finite temperature QCD and SQCD in the conformal window 
[27, 26, 25]. We have analyzed the dependence on p of thermodynamic quantities such as 
the entropy, energy and free energy density and of the jet quenching parameter and drag 
coefficient for a probe parton moving in the dual plasma. While in the "SQCD dual" 

q; a difference of [18] w.r.t. our lOd case is that we have total screening while, in the former paper, the 
screening only happens above a certain energy. 

13 It is not sensitive to the fact that the UV LST completion is a non-local theory. 
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the last two features are independent on p, in the "QCD dual" all the above quantities 
have been shown to be monotonically increasing functions of p. For the case of the jet 
quenching parameter, this seems to refine the possibility of using it as a direct measure of 
the temperature of the RHIC plasma, since it depends on the effective degrees of freedom 
(that, in turn, are also a function of T, albeit an unknown one). On the other hand, as 
a general feature, we have found that the dependence of all the observables on p is quite 
mild, so that the numerical results for the jet quenching parameter and the relaxation 
time are always very similar to the Af = 4 ones. This provides a first explanation of 
why the latter are so similar to the experimental results at RHIC. As far as we know, 
our results are the first example of unquenched effects of flavors in finite temperature 
holographic duals of gauge theories. 

As a second class of backgrounds, we have studied the lOd solution argued to be 
dual to a SQCD-like theory with Nf = 2N C flavors found in [29]. Its main features 
are that it exhibits total screening, the jet quenching parameter is zero and the drag 
coefficient is non vanishing. We have argued that these properties are common to all the 
finite temperature theories having LST as UV completion. This fact seems to force the 
interpretation that all these solutions in the UV give a description of the LST physics 
and so should not be reliable as realistic plasma duals in this regime. On the other 
hand, the expression of the drag coefficient turns out to be very similar to the M = 4 
one, underlying again the fact that dynamical quarks seem not to change drastically the 
physics of the plasmas. 

Both settings we have analyzed have some important problems. The non-critical one 
is not well under control, while the lOd one is probably thermodynamically unstable. 
Nevertheless, the results they provide are quite interesting and it should be worth in- 
vestigating other physical properties of finite p physics using these solutions. Of course, 
finding a finite temperature background which includes the backreaction of the flavors 
and does not share the problems of the above models would be the best starting point 
to study such properties. 
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A The black hole solution with spherical horizon 

In this section, we will introduce certain black hole metrics which generalize the solutions 
found in [29] for Nf = 2N C , Eq. (3.28). These black holes are characterized by an S* 3 
horizon and their temperature is always greater than the Hagedorn temperature of LST. 
The string frame metric of the compactified black hole solutions is 



ds = e 



- Tdt 2 + R 2 dQl + f Nc ^ ^ dr 2 + N c a' (-(d9 2 + sin 2 6d<p 2 ) + 

1 /„^Q .Q^„~Qn 1 



R 2 + N c a' 

+ -^—^(d9 2 + sm 2 dd{p 2 ) + + cosddip + cosddCp) 2 ^ , (A.l) 
with 

T = 1 - e 2ro ~ 2r , (A.2) 

so the horizon is placed at r = r . There is also a corresponding extremal solution which 
is recovered by setting T — 1. The £3 metric is written as 

dVt\ = dx\ + sin 2 x\dx\ + sin 2 x\ sin 2 X2 dx\ . 

The dilaton is linear 

$ = $ +r (A.3) 

and F( 3 ) is given by 

F( 3 ) = Nh sin 2 x\ sin xi dx\ A dx 2 A dx3 
N, 



, ^sin 9d9 Ad(f + sin 9 d6 A d(pj A (dtp + cos 9 dip + cos 9 dip) . (A.4) 
The radius R of the S3 horizon is determined by the flux through it 

(A.5) 



Since the flux Nh is quantized, the radius of the horizon takes only discrete values. By 
taking the formal limit R — > 00 of the above metric we find the solution of [29] . 
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The temperature and entropy density of the above black holes are 



1 / N c a' m / N c a' m , k , 

r(fl) = 2«V 1 + ^ = r »V 1 + ^ > T °- (A - 6) 



A 8 8e 4 *°iV r 4 „/ R 2 \ 3/ \ ., T 

R m ^T^ T ( R ) > ( A - 7 ) 



4G {10) C(4-0 V^ 2 + ^ c «' 
where To denotes the temperature of the black hole in [29]. Note that the temperature 
increases as the radius R of the horizon decreases. There is a hint that the above black 
holes may be thermodynamically unstable. Let us consider the above family of solutions 
in the limit R — > 00 and assume that they are a good approximation to black holes with 
a flat R 3 horizon like the solution in [29]. Then, as the black hole temperature increases 
as R decreases, (A. 6), we can see that the entropy s ~ R 3 is a decreasing function of the 
temperature, which is a sign of thermodynamical instability. 



A.l The qq potential 

In this section we calculate the potential between a probe quark and a probe antiquark 
in the above black hole backgrounds, as we have done in Section 3.1. The Wilson loop for 
two non-dynamical quarks separated by a distance L should be evaluated by the action 
of a fundamental string that is parameterized by t = r, x\ = a,r = r(a). Using the string 
frame metric 



ds 2 



- Tdt 2 + R 2 dVi\ + f Nc f t T- x dr 2 + N c a' (\{d6 2 + sin 2 8d<p 2 ) + 
' (dd 2 + sin 2 QdCp 2 ) + ^(dip + cos ddif + cos 6d(p) 2 j 



4-e 



(A.8) 



we find that the length L and the energy of this string upon subtraction of the infinite 
masses of the two quarks are given by 

L(r„) = 2^ r - = 2^Vl - e2ro-2r. K[e 2ro ~ 2r *], (A.9) 



HQ. 



(A.IO) 



1 VI - e 2r °~ 



•2r* ~2 



L^o 



2 2 Vl 3 ^ 2 



■ds - (1 



where r* is the minimum value of the radial coordinate reached by the string and 7 

R 2 N c a' 
R 2 +N c a> ■ 
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The length L(r*) is a monotonic function that goes from a minimum value of 
attained at r\ = r to the limit value of for r* — > oo. 14 The energy E(r*) vanishes 
for r* = ro <-> L = but is otherwise positive, which means that the quark and antiquark 
are always screened, as for the black holes (3.28). 

B Total screening in GTV black holes 

In this appendix, we will study the quark- antiquark potential in the black hole back- 
grounds found in [32]. We will find that, as in the case of the Nf = 2N C black hole, quark 
and antiquarks are always totally screened. Furthermore, this seems to be a generic 
feature of black hole backgrounds associated to systems of wrapped fivebranes. Indeed, 
one can repeat the analysis for the finite temperature solutions found in [44] and find the 
same result. As it transpires from the analysis below, this is related to the asymptotics 
of the solution both close to the horizon and far away from it. 
In the notation where the string frame metric reads 

during = -G 00 (s)dt 2 + G X] ^dx 2 + G ss (s)ds 2 + G XTXT {s)dx 2 T , (B.ll) 

define 

f 2 (s(x)) = Goo(s(x))G xn (s(x)), 

g 2 (s(x)) = Goo(s(x))G aa (s(x)). (B.12) 
The distance between quark and antiquark is given by the general formula [42, 43] 

2/" gww* (B.13) 



so 



In the case we are interested in, the string frame metric reads [32] 

ds 2 = e*< r > [ - v(r)dt 2 + dx 2 + u(r-y l dr 2 + R 2 (r)(d9 2 + sin 2 6d V 2 ) + l\ + e 2 + e 2 ] , 
so that 

f 2 (s(x)) = G 00 ( S (x))G :C||:C|| ( S ( a ;)) = ^(x))e 2 ^)), 

g 2 (s(x)) = G 00 (s(x))G s Mx)) = e 2 ^. (B.14) 

Close to the horizon, which we can set at r = [32] 

u = u' h r + 0(r 2 ), $ = $ ft + C(r), R = R h + 0(r), (B.15) 

14 For these black hole solutions with spherical horizon, x\ is an equatorial angle and L is then the 
angular separation between quark and antiquark along the equator. Note also that the limiting value of 
L is also less than tt. 
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whereas for r — > oo 

u=l + 0(e~ r ), $ = $ 00 + ^-ilogr + C(r- 2 ), R-V^- (B.16) 

Due to the asymptotic behavior of the above functions, the integral (B.13) is convergent 
in the limit s\ — > oo. 

We would like to show that L(s ) vanishes in the limit s — > 0. To this end, let us 
introduce the function 

T(s) = 1 , (B.17) 

y/u{s) (i/(s)e 2 *W-z/(s )e 2 *( s o)) 

that is useful to write L(sq) in a compact form, 

roc 

L(s )=2 \M^o) e* (so) T(s) rfs . (B.18) 

The function T(s) is analytic in [l,oo[ and its integral is convergent. Therefore, if we 
split the integral in (B.18) into two integrals from s to 1, and from 1 to oo, the latter 
simply vanishes in the limit s — > 0, as long as a/ v(sq) ~ y/so. Then, 



L(s ) = 2 I yM^)e^T(s)ds 

J So 



f 1 ds 

ryVso / --. ~ -v^o lQ g s o ■ (B.19) 

In summary, in the limit So — > 0, 

L(s ) ~ -^logso^O . (B.20) 

This behavior of the length functional is the same for all black holes with LST asymptotics 
[32, 44] and differs from other more familiar cases where L(s ) actually diverges in the 
same limit. It is also drastically different from the behavior in the zero temperature 
solutions [30, 46] where L(s ) actually diverges and there is linear confinement. 

We will now study the energy of the Wilson loop in the limit s — > and show that it 
is positive in a neighborhood of zero. The energy of the regularized Wilson loop, which 
means that we are subtracting the masses of the quark and anti-quark, is given by the 
following formula 



E(s ) = 2l I = ff= =-l\g(s)ds-2 / g(s) ds. (B.21) 



2 



so 



vm 2 - TOT 
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Given the asymptotics at infinity of the functions involved (B.16), the above integral is 
actually finite. Since the first integrand is positive 



E(s ) > E{s ) = 2 [ | - ^ = - 1 | g(s) ds-2 [° g{s) ds . (B.22) 

Thus, in the limit so 0, 

E(s ) = 2 r° ( - /(S) - 1 ) g(s) ds - 2 f ° g{s) 

4 Wf(s) 2 -f(*o) 2 J Jo 



ds 



2 / (u(s) T(s) e*^ - 1) e^^ ds-2 e^ s) ds 

Jso Jo 

2s ° e$(so) / 2 (t^t -1 ) rft " 2s ° e$(so) 

2v/2 + 2 log(l + V2) - 4) s e* (so) + O(sg) 



w 0.591 174 s e* (so) > . (B.23) 

The above is not enough to show that for this class of black holes quark anti-quark 
pairs are totally screened. In fact, it could be possible that L(s ) vanishes again and 
that around this value of so the energy is actually negative. However, numerical analysis 
shows that the length L is approximately a monotonic increasing function and does not 
vanish again. Furthermore, one can show that in the limit s — > 00 L(s ) goes to a 
constant value determined by the slope of the dilaton 

L(s )^^ <-> $(s)~Zs s^oo, (B.24) 

as we found in the Nf = 2N C black hole. Numerical analysis also confirms that the energy 
of the Wilson loops is always positive for s > 0. 

It would be interesting to check explicitly whether black hole solutions for Nf 7^ 2N C 
exhibit total screening as well. This fact seems to be determined by and intertwined 
with the asymptotics of the solution close to the horizon (B.15) and also the large r 
asymptotics (B.16). The fact that the energy of the Wilson loop is positive for small 
enough L and that the length presumably does not vanish again are strong hints that 
total screening might be a feature of general fivebrane or LST black holes. Indeed, one 
can perform the same analysis for the black holes found in [44] that have essentially the 
same asymptotic behavior as [32] and again find total screening. 
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